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Direct comparison is made of the steady-sates and coarsening dynamics in a local system and its 
nonlocal generalization. The example system is the surface of a solid film in a strong electric field; 
the morphological evolution of the surface is described, in the long-wavelength approximation, by 
the amplitude PDF for the film height function. It is shown that the amplitude of the steady-state 
and the coarsening rate of the surface structure are very sensitive to the radius of the long-range 
interaction, and that both quantities increase as the radius decreases. 
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Nonlocal pattern forming systems with long-range 
interactions, described by a nonlinear partial integro- 
differential equations, are ubiquitous in science and en- 

S ' reering. One can mention Rayleigh-Benard convection 
, magnetoconvection Q, surfactant-mediated interfa¬ 
cial flows Q, instabilities in plasma ii , flow of a film 
down an inclined plane in the electric field 0 , vibrated 
layers of a granular material or viscoelastic fluid @ , evap¬ 
oration of liquid films Q , and reaction-diffusion systems 

[Mini- 

In materials science, evolution of phases of a binary 
alloy is described by a Cahn-Hilliard (CH) equation for 
an order parameter; the impacts of long-range interac¬ 
tions on coarsening of the order parameter have been 
studied (l2l - [T5l| . Also, guided by their analysis of Asaro- 
Tiller-Grinfeld instab ility in heteroepitaxial solid films, 
Kassner and Misbah proposed the “generic” nonlo¬ 
cal amplitude equation, derived in the long-wave limit: 

ht = —ah + hxx + H [hx\ — -I- /3 {H [hx])^ , (1) 

where h{x) is the film height above the substrate, a, (3 
are parameters, and 

H[f[x)\ = 


is the Hilbert transform on the real line (p.v. stands for 
Cauchy principal value). Eq. ([1]), being relevant to sys¬ 
tems where the linear dispersion relation is quadratic, is 
complementary to the Kuramoto-Sivashinsky equation. 
Although these authors computed Eq. (P) and deter¬ 
mined that it describes a perpetual coarsening, the ac¬ 
count of their investigation is very brief (one figure dis¬ 
playing a qualitative behavior); in particular, they did 
not investigate how the non-local terms affect the coars¬ 
ening exponents. Nor were the effects of the finite radius 
of the long-range interactions studied. 

In this short note, we directly compare coarsening in 
the local and nonlocal systems described by a long-wave 
evolution equations for the film height. These equations 
stem from the consideration of the surface electromigra¬ 
tion [isl - f^ . an intrinsically nonlocal effect. Assuming 
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FIG. 1: (Color online.) Steady-state surface profiles from 
the evolution of the one-wavelength, small-amplitude sinu¬ 
soidal perturbation. This initial perturbation is shown by the 
dashed line, n = 2, 5, 200 curves are the steady-state profiles 
computed from Eq. (jd)). For comparison, the dash-dot line 
is the asinkmaxX curve, where a is the height of the n = 2 
curve above h — 0 level. The steady-states were checked by 
computing the evolution towards equilibrium of the perturbed 
profiles. 


the constant electric field Eq parallel to the substrate 
and using the local aproximation for the field on the film 
surface, E = Eg cos 0 (where 9 is the surface orientation 
angle), the local, conserved evolution equation reads [^ : 


-Bhxx. + A |-^ + M'(0) (^hx - y) } 

— Bhxxxx “ 1 “ AM (^O^hxx Ahxhxx 

- ^AM'mlhx., (3) 



where i? > 0 is the “Mullins number” characterizing the 
strength of the natural surface diffusion, A > 0 is the 
electric field strength parameter, and M'(0) < 0 is the 
derivative of the (anisotropic) adatom’s diffusional mobil¬ 
ity at the planar surface hx = 0. Under (dimensionless) 
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Eq. ([3|) a random, short-wavelength initial deformation 
of the infinite planar surface perpetually coarsens - the 
size L of the structures increases as a power law in time 
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kmax = ^ is the most dangerous wavenum¬ 

ber from the linear stability analysis of Eqs. and 
and Xmax is the corresponding wavelength. The 
boundary conditions at a: = 0, Xmax are periodic. R 
in Fig. [T] is represented by one half of the number of 
the grid points, n, over which the integral in Eq. 
is computed (that is, to compute the transform at Xi, 
the points a;_„, x_n+i,..., cci-i, x^+i,..., Xn_i, are 
used). R = oo (Eq. [2]) corresponds to all points in¬ 
volved, n = 200. It can be seen that when R = oo the 
steady-state amplitude is only slightly smaller than the 
one from the local Eq. ®, and the deviation from the 
local equation increases as the radius decreases. To re¬ 
affirm these results, in Fig. [2] we plot the steady-state 
amplitude and maximum slope vs. R; both quantities 
decrease fast initially and level off already at n ^ 75. 
And, from the monotonicity of the graph of the ampli¬ 
tude of the steady-state vs. the wavelength (Fig. [S]) one 
concludes that coarsening is uninterrupted in a nonlocal 
system 


FIG. 2: The amplitude and maximum slope of the steady- 
state profile vs. (twice) the interaction radius. 

Guided by Ref. [l^ , the nonlocal generalization of this 
equation is: 
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R is the interaction radius. We corupute Hilbert trans¬ 
form using fast method from Ref. [l7[. The reason we 
focus on Eq. ®, rather than 
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is because the latter equation describes the unlimited 
slope growth, and thus it does not allow to quantify the 
long-time effects of the nonlocal terms. The unlimited 
slope growth is inconsistent with the fully nonlocal com¬ 
putations of the morphology evolution using the moving 
boundary problem (where the electric field is com¬ 
puted in the bulk of the film using the Laplace equation 
for the electrical potential). 

Fig. [T] shows the steady-state profiles from the initial 
condition /i(a:, 0) = 0.01 sinkmaxX, 0 < x < Xmax, where 



FIG. 3: Amplitude of the steady-state surface profile vs. its 
imposed wavelength A. n = 5. 

Next, we performed computations of coarsening using 
Eqs. ([3]) and (H)). All such computations are done on 
the domain 0 < x < 20Xmax with the periodic boundary 
conditions. The results are averaged over five runs with 
a different random initial condition, and numerical con¬ 
vergence is checked on large grids. Here R = oo (Eq. [2|) 
corresponds to n = 400. In Figures SK a,b) L is the hor¬ 
izontal scale of the surface structure, e.g. the mean size 
of a hill at its base. One can see that coarsening is very 
slow for large R (or n) and speeds up when R decreases, 
with the rate approaching one from the local equation 
as i? —0. The major speed-up occurs in a quenching 
fashion when n decreases from 8 to 5, and from there 
the speed-up is gradual. Also the coarsening rate is very 
weakly sensitive to n for n >ss 10 (Fig. |D(b)). 
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The coarsening exponent in the local model, Eq. 
initially is of the order reported in Ref. [l^ for a similar 
local model (0.4 vs. 0.3), then it decreases sharply for 
the rest of the evolution. In the nonlocal model, as R in¬ 
creases, the initial regime is confined to the progressively 
shorter time intervals, and then coarsening becomes log¬ 
arithmically slow. The initial power law coarsening is 
matched to the logarithmic law by a power law with a 
smaller exponent, as shown in the inset of Fig. UKa). 
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FIG. 4: (Color online.) (a), (b): Length scale of the coarsen¬ 
ing surface structure vs. the time, (b) shows the zoom into 
(a). The thin black lines in (a) are the power law fits to the 
data. 
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